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In this paper we examine the design of scalar and 
multi-variable feedback control systems for the GE 
T700 turboshaft engine coupled to a helicopter 
rotor system. A series of linearized models are 
presented and analyzed. Robustness and perform- 
ance specifications are posed in the frequency 
domain. The LQG/LTR methodology is used to 
obtain a sequence of three feedback designs. Even 
in the single-input single-output case, comparison 
of the current control system with that derived 
from the LQG/LTR approach shows significant 
performance improvement. The multi-variable 
designs, evaluated using linear and nonlinear sim- 
ulations, show even more potential for performance 
improvement. 


1. INTRODUCTION 

In this paper we summarize, til, three distinct 
feasibility studies related to the design of feedback 
control systems for a model of the GE T700 
turboshaft engine coupled to a helicopter rotor 
system. The present control system on the T700 
engine uses a single input the fuel, and was 
designed using classical single-input single-output 
(SISO) techniques. We explore the potential 
advantages of using more sophisticated com- 
pensators, derived using the Linear-Quadratic- 
Gaussian with Loop-Transfer-Recovery (LQG/LTR) 
design methodology, both in the SISO case and in 
the multiple-input multiple-output (MIMO) case. 
In the MIMO case we use the dynamic coordination 
of both fuel and variable compressor geometry to 
control two outputs of interest. 
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To the best of our knowledge this is the first study 
dealing with the application of multi-variable 
design concepts to a turboshaft engine. On the 
other hand, the modern multivariable control of 
turbofan engines has received a great deal of 
attention. The book by Sain et al, [2], contains a 
variety of design studies on the F-100 turbofan 
engine; other pertinent references are [3] to [1 1]. In 
particular, feasibility studies using the LQG/LTR 
design methodology have been reported for the F- 
100 engine in [3] and [4], the GE-21 engine in [5], 
and the GE-16 engine in [6]. There seems to exist 
widespread agreement that the dynamic coordin- 
ation of fuel with several engine geometry vari- 
ables will result in future multi-variable feedback 
designs that will improve engine efficiency, result 
in more rapid thrust response, tighter control of key 
temperatures and pressures, and improved stall 
margins. 

The dynamic models used in this study include the 
interaction between the turboshaft engine and the 
helicopter main-rotor and tail-rotor dynamics. As 
explained in Section 2, we included the engine- 
rotor dynamic interactions in our model because 
the bandwidth specifications, that we have imposed 
to carry out our feasibility studies, were larger than 
those of the production design, and consequently 
the resonances associated with the main and tail 
rotor dynamics had to be included in our model. Dn 
the other hand, precise knowledge of such re- 
sonances is not available. For this reason, we have 
estimated engine-rotor model errors in the 
frequency . domain, and imposed stability-robust- 
ness specifications, 90 as to account for such 
modeling errors. We do not claim that we have 
captured all relevant high frequency modeling 
errors;, nonetheless, a similar stability-robustness 
analysis will have to be carried out in a more 
realistic application. 

We present evaluations of three distinct feasibility 
studies for the engine-rotor system. Design A is a 
SISO design using the LQG/LTR method. In Design 
A we use only the fuel to control the free (power) 
turbine speed. We compare the "sophisticated" 
r Design A with the existing production design, and 
'^demonstrate improved performance. Thus, there 
exists potential performance payoff in using, even in 
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a SISO setting, dynamic models of greater fidelity 
and more sophisticated compensator designs. 

Design B is a MIMO design, and it is used to 
demonstrate the advantages of usinjr an additional 
control variable. We use the dynamic coordination 
of both fuel and variable compressor geometry to 
independently control the free turbine speed and 
the gas generator speed. We compare the MIMO 
Design B to the SISO Design A with respect to their 
disturbance rejection properties. We show that 
Design B is superior in the sense that dynamic 
modulation of the variable geometry control is 
utilized to reject disturbances. The same disturb- 
ance in the SISO design would have to modulate 
the fuel, thereby possibly decreasing the engine 
fuel efficiency. ' 

Design C is a different MIMO design. Like Design 
B we again use both the fuel and variable geometry 
as dynamic controls. However, in Design C the 
outputs that we wish to control are the free turbine 
speed and the inter-turbine gas temperature. 
Precise control of temperature is necessary when 
the engine operates at nigh-power conditions so as 
to prevent damage. With this choice of controls and 
outputs the engine-rotor open-loop dynamics 
exhibit a non-minimum phase zero at about 0.2 
rad/sec. The presence of this non-minimum phase 
zero imposes limitations with respect to the 
command-following and disturbance-rejection per- 
formance of the feedback system. Nonetheless, we 
demonstrate that slow temperature trim 
commands, useful for dynamically improving 
engine efficiency, can be reasonably followed. 

At this point it is important to stress that the 
results presented in this paper only represent 
feasibility studies, and more work is needed before 
the LQG/LTR based compensators are implemented 
in a working control system. Although we have 
used non-linear simulations to evaluate the 
designs, we did not test them over the full envelope 
of possible operating conditions. It is likely that 
gain-scheduling will have to be used to develop 
designs that maintain improved performance and. 
stability over the full operating envelope. In spite 
of these limitations, the results demonstrate 
significant advantages of using multi-variable 
control for turboshaft engine applications, and 
illustrate the types of compensators that would 
result from the LQG/LTR design methodology [12] 
to [17]. 

The remainder of the paper is organized as follows. 
In Section 2 we present a discussion of the GE T700 
turboshaft engine dynamics and its dynamic 
coupling to the nelicopter rotor system, including a 
discussion of the nature of the linearized dynamics. 
In Section 3 we present an analysis of the linearized 
dynamics for tne three designs, in terms of poles 
and zeros, and frequency domain singular value 
plots. In Section 3 we also quantify the modeling 
errors in the rotor dynamics so that we can im- 


pose stability-robustness specifications in the fre- 
quency domain; in addition, we summarize the 
idealized performance specifications that we im- 
posed for our feasibility studies. In Section 4 we 
first present a brief overview of the LQG/LTR 
design methodology which was used for deriving 
the Designs A, B, and C. In Section 5 we sum- 
marize the characteristics of all three designs in the 
frequency domain by presenting the shapes of the 
singular values of the loop, sensitivity, and closed- 
loop transfer function matrices vs. frequency. 
Then, we evaluate the transient performance char- 
acteristics via simulation. Section 6 summarizes 
the conclusions. The appendix contains the state 
equations and the numerical values of the open- 
loop dynamics in terms of the A, B, C, D matrices of 
the state-space models. 


2. SYSTEM DESCRIPTION AND MODEL 
FORMULATION 


2.1 System Description 

A conventional helicopter, as shown in Figure 1, 
utilizes a single main rotor, primarily for lift, and a 
tail rotor for torque reaction and directional control 
in the yaw degree of freedom. The main and tail 
rotor systems are directly coupled to two turboshaft 
engines through gear reduction sets and shafting! 

The main and tail rotor systems are composed of in- 
dividual blades which are simply airfoils that 
royide lift and/or thrust. The pilot maneuvers the 
elicopter by modulating the available lift/thrust 
from the rotor systems. A maneuvering demand 
from the pilot is equivalent to producing a load dis- 
• turbance on the rotor systems. Load disturbances 
may also emanate from other sources such as wind 
gusts. The incorporation of a "fast" or "tight" 
engine speed control capable of rejecting rotor 
system load disturbances will be reflected in 
increased helicopter maneuvering capability. This 
increased maneuvering capability must be accom- 
plished, however, without exciting coupled 
engine/rotor system complaint dynamics that are 
present. 

The turboshaft engine utilized in this study is the 
GE T700 engine, as representative of a recent 
technology engine in current production. A 
simplified cross-section of the GE T700 is shown in 
Figure 2. The gas generator sustains the gas tur- 
bine cycle, while the free turbine performs the role 
of extracting energy. It is the free turbine, when 
.directly coupled to the helicopter rotor system, that 
recovers the useful work of the gas turbine cycle. 
The responsibility of the gas generator is to provide 
the power demanded by the helicopter rotor sys- 
tems at a specified free turbine speed. The turbo- 
shaft engine control system must insure that the 
power demanded by the helicopter rotor system is 
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supplied by the engine while simultaneously 
insuring that the engine operates efficiently over a 
wide range of ambient conditions and prevents 
destructive stall phenomenon, turbine overtem- 
peratures, overspeeds and excessive shaft torque. 
On the GE T700 two control variables, fuel flow 
and compressor variable geometry, can be utilized 
to meet these objectives. In the current control 
system only the fuel is controlled in a closed-loop 
sense. The compressor variable geometry is sched- 
uled in an open-loop manner. 

The operational requirements for the turboshaft 
engine are summarized below: 

The engine must: 

1. Maintain constant free turbine speed in the 
presence of load disturbances to the helicopter 
rotor systems, 

2. Not provide input energy to excite coupled 
engine/drive tram resonant modes, 

3. Maintain adequate stall margin, 

4. Limit turbine inlet temperatures, speed and 
torque, and 

5. Operate at peak efficiency. 

The remainder of this paper overviews the results 
of research 11 ] undertaken to examine the 
feasibility of the LQG/LTR control design 
methodology in achieving the above goals. 

2.2 Model Formulation 

Turbine engine dynamics are described by complex 
non-linear equations relating state variables x(t), 
control variables u(t), output variables y(t) and 
ambient variables in the form 


d/dt x(t) = f (x(t), u(t), 0 ) 

(la) 

y(t) = g(x(t),u(t), 0 ) 

(lb) 


The linear, time-invariant, constant coefficient 
model utilized in the subject research is of the form 


d/dt 8 x(t) = A 8 x(t) + B 8 u(t) (2a) 

Stft) = C 8 x(t) + DSu(t) ( 2 b) 


where 
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The objective in model formulation is to establish a 
nominal representation of the open-loop system, or 
plant, in the low frequency region where per- 
formance specifications are imposed. To achieve a 
practical and implementable design, the coupled 
engine/helicopter dynamics in the 0-40 raa/sec 
frequency range were examined for inclusion in the 
nominal linear model. The two available control 
variables, fuel flow, Wf, and compressor variable 
geometry, V g , were included in tne model repre- 
sentation to provide independent control of two 
output variables (to be discussed later). 


The low frequency (< 10 rad/sec) GE T700 engine 
dynamics are dominated by the gas generator and 
free turbine dynamics. Pressure and temperature 
dynamics appearing in the flow equations are 
typically "fast" for a small turboshaft engine and 
are included in the model only as outputs, thus 
neglecting their dynamics. Inter-turbine gas tem- 
perature, T 45 , was included in the model as an out- 
put as it is often desired to control that variable. 
The reduced engine state vector for design purposes 
is thus given simply by the two turbine speeds: the 
|as generator speed, N g , and the free turbine speed. 


The state variables are associated with energy 
storage elements and are temperatures, pressures 
and inertia terms for a gas turbine system. The 
control inputs are fuel flow and variable 
geometries. The outputs can be turbine speeds, 
pressure ratios and gas temperatures. The ambient 
variables are ambient pressure and temperature 
ratios. 

Prior to formulation of the control problem, the 
non-linear dynamic description must be converted 
to a linear dynamic model pertinent to operation 
about an equilibrium operating condition. The 
equilibrium condition is characterized by 0 and the 
steady-state values of the state, control and output 
variables (xo, Uq, yJ. 


The helicopter drive train compliant dynamics 
must be represented in the system model because 
they are present within the engine response 
bandpass. A representative helicopter drive train 
is shown isometrically in Figure 3. A simplified, 
lumped parameter, spring-mass-damper represen- 
tation of the system is shown in Figure 4. 

The turboshaft engines are coupled to the 
-helicopter drive-train model as shown in Figure 5, 
which is a block diagram representation of the 
coupled system. Note that the only coupling is 
through Q p , which is the gas torgue generated by 
the gas generator and applied at the power 
turbines. The state variable representation of 
Figure 5 is given in Appendix A. The variable def- 
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initions for Figure 5 with units are included in 
Table 1. The full model also includes T4.5 sensor 
and Vg actuator dynamics. 

A series of linear models were generated re- 
presenting a range of operating conditions. The 
operating conditions are defined by power level or 
percent of design gas generator speed. 


3. SYSTEM DEFINITIONS, ROBUSTNESS 
REQUIREMENTS, AND DESIGN 
SPECIFICATIONS 


3.1 System Definitions 

Three distinct open-loop system definitions are pre- 
sented for utilization in controller design. The first 
(System A) is a conventional SISO system with 
scalar fuel flow control (Wf) and a single scalar 
output, the free turbine speed, No- With two 
available control inputs, fuel flow, Wf, and variable 
geometry position, V g , control over two distinct 
output variables is realizeable. Accordingly, the 
second system definition (System B) explores the 
simultaneous control of free turbine, speed, N p and 
gas generator speed, N g . The third system 
definition (System C) represents an exploration 
into the simultaneous control of free turbine speed 
and inter-turbine gas temperature, T4.5. 

The control of power turbine speed is required to 
satisfy the fundamental system requirement of a 
commanded power supply to the helicopter rotor 
systems. In System B, the simultaneous control of 
the power turbine and gas generator speeds was 
undertaken to explore the utilization of this control 
system definition for both input and output 
disturbance rejection as compared to the SISO 
controller. The simultaneous control of turbine 
temperature in System C allows a potential handle 
on dynamic engine operational efficiency and 
provides some latitude in temperature limiting. 

The input-output definitions are summarized in 
Table 2 along with the operating conditions for the 
linear dynamic model utilized for each design as 
denoted by the power level or % N g . The 90% N g 
design model was chosen for the SISO design 
(System A) because it is representative of normal 
operating power. The 83% N g design model level 
was utilized for MIMO System B to examine the 
implications of a MIMO control law at a low power 
level. The 90% N g design model was chosen for 
MIMO System C because it is, as mentioned pre- 
viously, representative of normal operating power. 
The linearized equations for each system are 
summarized in Appendix B. 

The poles and zeros of the design models for each 
system definition are tabulated in Table 2. The 


zeros for the multi-variable system definition are 
transmission zeros [13]. Note that System C has a 
non-minimum phase zero at 1199 rad/sec which pre- 
sents a generic performance limitation [20] that 
will become clearly evident in the controller design 
results section. 

Figures 6 to 8 display the singular values of the 
three open-loop Systems A, B, and C respectively. 
Thus, we plot the singular values [17] 


Oi [Gp (jci»)l (3) 

where 

G p (s) = C (s I - A)-i B + D (4a) 

X(s) = G p (s)u(s) (4b) 

is obtained from the linear models in Appendix B. 
All three frequency plots show the effects of main 
rotor and tail rotor resonances. Also note that (see 
Figure 8) the presence of the non-minimum phase 
zero for System C manifests itself in a small 
minimum singular value at low frequencies, as 
compared to Systems A and B. 

3.2 Robustness Requirements 

The dominant high-frequency uncertainty in the 
linear model is in the description of the helicopter 
rotor dynamics. The rotor system lumped para- 
meter model does not portray the functional 
relationships of the main rotor spring and damping 
coefficients with helicopter flight condition, rotor 
coning angle, etc. The posing of a maximum 
realizeable range of coefficient variation, while not 
capturing the explicit functional relationships, will 
acknowledge their presence and provide the basis 
for a conservative, stable design. 

Consider the actual plant, G (s) to be related to the 
linear model representation, by the expression 

5 (s) = L (s) G p (s) (5) 

Equation (5) relates the model uncertainty 
quantified by L(s) to the system output variables 
as shown in Figure 9. 

For the error defined by Equation (5), an output 
feedback system is guaranteed to be stable if the 
inequality 

Oma* t LCI«)-I ] £ o roin [ I -I- (Gp(jo>) K(jo>)) -1 J (6) 

is satisfied for all o> [18, 19], where K(jo>) is the 
compensator transfer function matrix. The error 
matrix, L0‘o), will be assumed to be of the form 

L(jw) = I + E(jo>) = I + diag [e r (o>), e e (co)] (7) 
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where 

e r (o) = rotor system error 

e 4 (o) = engine system error. 

The above error structure reflects the uncertainty 
of the high frequency rotor system description to 
the power turbine speed output. The error in the 
low frequency engine system will be assumed small 
and we will let e e (co) = 0. 

The construction of e c (co) to quantify rotor system 
parametric variations i9 best visualized on a polar 

( dot (Nyquist Diagram) of the rotor system open- 
oop transfer function as shown in Figure 10. The 
rotor system dynamics are present in the system 
transfer function matrix for free turbine speed 
output and for both control variables Wf and V g . 
Realizeable rotor system variations perturb the 
nominal system representation. A circle of radius r 
and center coincident with the nominal model 
encompassing the family of perturbed plants at 
several frequencies is also shown in Figure 10. The 
function relationship of the magnitude of r with 
frequency establishes e r (o)). The function e r (u) 
thus quantifies the realizable magnitude 
variations of the rotor system representation so 
that a compensated system can be designed that 
does not realize a change in the number of 
encirclements of the critical point on the Nyquist 
Plot. A change in the number of encirclements is 
indicative of instability, and will be avoided if the 
inequality presented by Equation (6) is satisfied. A 
logarithmic plot of e r (o>) is shown explicitly in 
Figure 11. 

3.3 Design Specifications 

The implication of specifications is to achieve good 
performance in terms of 

1. command following, 

2. disturbance rejection, and 

3. insensitivity to modeling error through 
the introduction of feedback. 

For a feedback system, as shown in Figure 12, the 
maximum error at a given frequency, o> 9 , for unit 
magnitude commands and output disturbances is 
given by 

|et J =l/o min [I + T(jo) 0 )] (8) 

where 

T(ju 0 )= G(ju 0 ) K(j« 0 )- 

If OminXO^o) ► 1. then Equation (8) is nearly 

Ulj = l/Omml(jWo)- (9) 


Intuitively, by making o m ; n T (jco 0 ) "large" over a 
wide frequency range, we can both reject output 
disturbances and follow commands with small 
errors. This performance consideration must be 
tempered with a bandwidth limitation, (i.e. when 
°min X (jo>o) = 1) so that unmodeled dynamics do not 
cause instability. A desirable crossover frequency 
for this system, as derived from pilot evaluations 
[21] is about 10 rad/sec. To provide maximum 
command following and disturbance rejection, it is 
desired that 

Omin X (j«) > 20db V to ^ 1 rad/sec. (10) 

It is also required that all output variables have 
zero steady-state error to constant reference inputs, 
thus dictating integral augmentation. 

Figure 13 summarizes the frequency domain 
performance specifications. Robustness will be 
achieved through satisfaction of the inequality 
presented in Equation (6). 

The above specification should be viewed as 
tentative. It is well recognized by now that the 
presence of low-frequency non-minimum phase 
zeros (as we have for System C) represents a generic 
limitation in performance independent of the 
design methodology emjployed. As we shall see in 
the next section, we will not be able to meet the 
specification above for System C. 

The following sections present the design 
methodology and the controllers designed to meet 
the specifications. 


4. LQG/LTR DESIGN METHODOLOGY 
OVERVIEW 

The control structure to be utilized, including 
integral augmentation, is shown in Figure 14. The 
LQG/LTR compensation to be designed is given by 
Kn(s). The overall compensator, which includes 
integral augmentation, is defined by 

K(s) = (I/s) Kd(s). (11) 

For design purposes, the integral augmentation is 
considered part of the plant, thus we define: 

G a (s) = G p (s)(I/s). (12) 

The LQG/LTR procedure begins with the state 
description of the augmented plant given by 

d/dt X ft (t) = Aa Xa(t) + Ba Ua(t) (13a) 

Xa(0 = CaSaU) (13b) 

and thus 


Sa(s) = Ca(sI-AaH!a. 


(13c) 





The step-by-step LQG/LTR Cl 6, 17] design proced- 
ure to define the compensator Kd(jco) is: 

Step 1 : Shape the Kalman Filter Loop Transfer 
Function GkfCM given by 

GkfGu) = Ca (i«I* Aa) - 1 Ha (14) 

where Ha = Kalman Filter Gain Matrix using the 
free parameters L and p available in the Kalman 
Filter algebraic RTccatti equation 

Aa£ + E Aa + L LT-(l/p) E Cl Ca E = 0 (15) 

to yield the gain matrix 

Ha = (1/p) E Qa (16) 

so that the performance specifications posed in 
Section 3 are met. 

Step 2 : Solve the following algebraic Riccatti 

equation. 

KAa + AjK + qCjCa - KgaBjK = 0 (17) 

for q -* <» (sufficiently large) to yield the control 
gain matrix 

Ga = BjK. (18) 

As q oe> the LQG/LTR method guarantees that in 
the absence of non-minimum phase zeros 

Ta(s) = Ga(s) Kd(s) - Gkf (s). (19) 

Thus if we design GkfO'w) to meet the posed 
frequency domain performance specifications, and 
there are no non-minimum phase zeros in the 
system, then we can design a compensator, defined 
by 

Kd(s) = Ga(sL- Aa + BaGa + HaSa) “ l H*(20) 

by utilizing the asymtotic adjustment procedure 
(i.e. "cheap” LQG control problem) defined in Step 
2. The presence of a non-minimum phase zero 
within the desired bandwidth of tne system 
presents a generic performance limitation that 
cannot be considered an indictment of this 
methodology. The restrictions presented by a non- 
minimum phase zero within the desired system 
bandwidth will be demonstrated in the following 
section. 


6. CONTROLLER DESIGN RESULTS 


5.1 S1S0 Design A 

For the SISO system definition with scalar fuel 
flow control over the free turbine speed, the 


application of a generalization of the design 
methodology presented in Section 4 to the SISO 
case results in the loop transfer function shown in 
Figure 15. The system is shown to achieve the 
desired performance specifications. The com- 
pensator transfer function is shown in Figure 16. It 
is readily seen by comparison of the open-loop SISO 
transfer function of the plant shown in Figure 6 
with the compensator transfer function of Figure 
16, that the LQG/LTR design methodology 
perforins an approximate inversion of the open-loop 
plant [17]. Robustness is achieved as shown in 
Figure 17 for the modeling errors quantified in 
Section 3. 

A comparison using nonlinear dynamics was 
performed to determine the increase in system 
performance achievable with a LQG/LTR controller 
vs. a lower bandwidth, conventional controller. 
Figure 18 displays a typical helicopter transient, a 
30% load demand performed in 1 sec., for both a 
LQG/LTR and the current conventional controller. 
The fact that the LQG/LTR controller provides 
"tighter" power turbine speed governing is readily 
observed by comparing the power turbine speed 
deviations shown in Figures 18a and 18b. The 
improvement can be quantified by noting that the 
sensitivity transfer function, given by 

1/(1 + g(s) k(s)) (21) 

and shown in Figure 18c, which is indicative of the 
systems response to load disturbances, provides load 
disturbance attenuation over a much wider range of 
frequencies than does the conventional controller. 

5.2 MIMO System B Design 

The first MIMO design, providing the coordinated 
control of the power turbine and gas generator 
speeds using the fuel flow and variable geometry 
inputs, is presented to demonstrate not only the 
LQG/LTR design methodology, but to demonstrate 
how the coordinated control of several variables can 
provide a performance not realizeable with conven- 
tional scalar controls. 

The loop transfer function designed for this system 
definition is shown in Figure 19. Note that the 
performance specifications are met. Robustness, for 
the modeling errors quantified in Section 3, is 
achieved as shown in Figure 20. The MIMO 
compensator transfer function KqOu) is shown in 
Figure 21. Note, as in the SISO case, that the 
compensator, Kd(j&}), * 9 an approximate inversion 
of the open-loop plant shown in Figure 7. The 
closed-loop and sensitivity singular value plots are 
shown in Figure 22 and 23, . respectively. The 
closed-loop and sensitivity singular value plots 
demonstrate that good command following and 
disturbance rejection is realized. 

It is instructive to perform a comparison between 
this MIMO system definition and the SISO system 




definition in terms of the ability ofthese systems to 
reject disturbances. The basis for comparison will 
be to determine if it is possible, by using the co- 
ordinated control of several variables, to provide a 
svstem that exhibits disturbance rejection cap- 
ability without the extensive use of control energy. 
If we consider that the variable geometry input 
variable is available at no cost to the user, then any 
use of the variable geometry input can be consid- 
ered a savings of fuel. 

It is possible to incur a disturbance in the gas 
generator speed due to power extraction or engine 
inlet distortion. A linear simulation of the 
response of MIMO Design B to a step disturbance 
on gas generator speed is shown in Figure 24a. A 
linear simulation of the same disturbance to the 
SISO design A is shown in Figure 24b. Note that 
while both the MIMO and SISO systems reject the 
gas generator speed disturbance, the MIMO system 
rejects the disturbance rapidly and with no steady- 
state fuel cost as is incurred in the SISO system. 

The comparison between the MIMO and SISO case 
is instructive in pointing out that extended system 
performance capabilities are possible through. the 
coordinated control of several variables, provided 
that no generic limitations induced by non- 
minimum phase zeros exist. 


6.3 MIMO Svstem C Desig 


MIMO System C provides for the control of both the 
power turbine and inter-turbine gas temperature, 
using the fuel flow and variable geometry controls. 
A non-minimum phase zero is present at .199 
rad/sec in this design model due to the interaction 
of the variable geometry and airflow/temperature 
dynamics in the engine. The singular value loop 
transfer function for this design is shown in Figure 
25. Note that the posed performance specifications 
cannot be met. This is due to the presence of the 
non-minimum phase zero at .199 rad/sec which 
limits the frequency range for which Omj n T(jw) can 
be made "large”. The system is robust however, for 
the modeling errors defined in Section 3 as shown 
in Figure 26. The closed-loop and sensitivity plots 
are shown in Figures 27 and 28, respectively. Note 
that the effect of the non-minimum phase zero is 
demonstrated in all singular value plots. 

A trim signal on the inter-turbine gas temperature 
reference, a 10-second ramp of 20'T, as performed 
on a non-linear simulation is shown in Figure 29. 
The error magnitude between the temperature 
reference and the sensed temperature is deter- 
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mined by o m i n T(jo>) as presented by Equation (9). 
The 10-second ramp trim signal was chosen to 
represent the fact that the command error mag- 
nitude could be maintained at a relatively small 
level if the trim command signal is in the frequency 
range for which o m in T(jco) is ’’large". 


6.4 Design Summary 


Three important issues were demonstrated in this 
section. The first is that the LQG/LTR design 
methodology provided a systematic approach to 
frequency domain 'loop shaping’ for both the SISO 
and MIMO case. Additionally it was shown that 
the coordinated control of several variables can be 
utilized to provide performance not achievable with 
conventional scalar controls. Finally, the generic 
performance limitations of non-minimum phase 
zeros were demonstrated. 


6. CONCLUSIONS 


In this paper we have examined both SISO and 
MIMO designs for the control of a model of the GE 
T700 turboshaft engine including its dynamic 
coupling to the rotors of a helicopter. All control 
system designs were carried out using the so-called 
LQG/LTR design methodology. The results in- 
dicate that there is potential for significant payoff 
in command-following and disturbance-rejection 
performance, if realistic models of the engine-rotor 
system are used in the design process. It was also 
demonstrated that the dynamic coordination of 
both fuel and variable geometry controls results in 
superior performance. There was no particular 
difficulty in applying the LQG/LTR procedure to 
these designs, even when the plant had non- 
minimum phase zeros (Design C). In the latter 
case, the results were predictable and consistent 
with the limitations in performance inherent in 
non-minimum phase systems. 

We reiterate that our results should be viewed as 
feasibility studies. Much more work is needed to 
design a full envelope control system for the engine- 
rotor dynamic system. 
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TABLE 1. VARIABLE DEFINITIONS 


Dam m Main Rotor Aerodynamic Damping, fV-lba/RPM 
Dat ° Tail Rotor Aerodynamic Damping, ft.-Ibs/RPM 
Djifl <■ Main Rotor Damping, fL-lbs/RPM 
Jur a Main Rotor Inertia, It.-lb.-aee./RPM 
J T ■ Lumped Drive Train Inertia, ft.-lb.-see./RPM 
J g *> Cas Generator Inertia, ft.-lb.-sec./RPM 
Kiah * Main Rotor Spring Conatant, fL-lbs/RPM-see 
K th ** Tail Rotor Spring Constant, fL-lbs/RPM-aec 
N|u ° Main Rotor Speed, RPM 
Ntp * Tall Rotor Speed, RPM 
N, “ Gaa Generator Speed, RPM 
N p *■ Free Turbine Speed, RPM 
Qmr m Main Rotor Torque, ft-lbs. 

Q th ■ Tali Rotor Torque, R--lba. 

Q, «■ Cas Generator Cas Torque, fL-lba. 

Qp » Power Turbine Gae Torque, ft.- lbs. 

T 4.5 b Interturbine GasTemperature, degrees R 
T tc b Thermocouple Lag, 1/sec. 

Tvo 0 Variable Geometry Actuator Lag, l/aec. 

V ( b Variable Geometry Position, degrees 

V r » Variable Geometry InputCommand, degrees 

W f - Fuel Flow, IbThr. 


TAIL ROTOR 



Figure 1. Conventional Single Main Ilotor 
Helicopter 


’FUEL FLOW INPUT (CONTROL) 
VARIABLE / T Q COMBUSTOR 

6E0METRT (CONTROL) ' l - uw,ua u 



GAS GENERATOR 


Figure 2. GE T700 Simplined Cross Section 


TABLE 2. OPEN- LOOP POLE-ZERO 
STRUCTURE OF DESIGN MODELS 


SYSTEM 

DEFINITIONS 

POLES 

ZEROS 

SISO SYSTEM A 

OPERATING- 

626 

-7.37 

CONDITION: 

-3.64 

-.86 + ,6.84 

90% N g 

-5.07 + ,15.8 

-.18 + ,34.2 

CONTROL: W, 
OUTPUT: N p 

-.62 +,404 


MIMO SYSTEM B 

OPERATING 

-.482 


CONDITION: 

-1.92 

-.84 + ,6.85 

83% N g 

-10.0 

CONTROLS: W,.V g 

-5.07 + ,16.8 

-.18 + ,34.2 

OUTPUTS: Np.T 4 ,j 

- .62 + ,40.4 

MIMO SYSTEM C 

OPERATING- 

.626 

.199 

CONDITION: 

-2.2 

-.85 + ,6.85 

90% N g 

-10.0 

CONTROLS: W f ,V g 

-6.07 + ,15.8 

-.18 + ,34.2 

OUTPUTS: N p ,T 4 .j 

- .62 + ,40.4 




Figure 3. Helicopter nutor System Drive- 
Train Isometric Diagram 



Figure 4. Lumped Parameter Representation of 
Helicopter Ilotor System Drive Train 
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Figure S. Block Diagram Representation of Coupled 
Turboshaft Engine Rotor System 


(A) 



Figure 0. Magnitude Dode Plot of the Open-Loop 
Transfer Function for the SISO 
System A 


(A) 



Figure 7. Singular Value Plot of the Open-Loop 
Transfer Function for the MlMO 
System D 



Figure 8. Stngulnr Value Plot of the Open-Loop 
Transfer Function for the MIMO 
System C 


o Is) 


HAST 
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Figure 9. Modeling Error Definition 






la 



Figure 10. Depiction of Rotor System Unstructured 
Modeling Error 


(A) 



Figure 11. Bound on Multiplicative Error Versus 
Frequency 
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Figure 13. Frequency Domain Design Specifications 



A * -1 

Figure 14. Control Structure Utilized in Design Process 


(db) 



Figure 12. Feedback System Definition 


Figure 15. Magnitude of SISO Loop Transfer Function 
g(jcj) k(jti>) from LQG/LTR Design A 
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Figure 18. SISOLQG/LTR Compensator k D 0‘«J 
Magnitude vs. Frequency for Design A 


<db> 



(rsd/sec) 
LOG FREQUENCY 


Figure 17. Stability Robustness Check for Design A 



Figure 18a. Non-linear Simulation of Current Controller 
Compensated System Response to 30% 

Load Demand in 1 Second 



Figure 18b. Non-linear Simulation of SISO LQG/LTR 
Compensated System A Response to 30% 
Load Demand In 1 Second 





u 



Figure 18c. Comparison of the Sensitivity Transfer 
Functions of the Current and LQT/LTR 
Control System Oeslgn A 
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Figure 21. Singular Values of LQT/LTR Compensator 
K d (Jcj) for Oeslgn B 



Figure 19. Singular Values of Loop Transfer Function 
G(Ju) K(Jw) for Design B Using 
LQG/LTR Compensator 



Figure 20. Stability Robustness Check for Oeslgn B 
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Figure 22. Singular Values of Closed-Loop Transfer 
Function G(ju) K(jcj) II + G(Jo>) KOwJl -1 
of Design B 


o 



LOG FREQUENCY 


Figure 23. Singular Values of the Sensitivity 

Transfer Function (I + £<Jw) K(J<i>))-* 
of Oeslgn B 











(rtd/KC) 
LOG FREQUENCY 


Figure 20. Singular Values of Loop Transfer Function 
G(ju) K(Ju») for Design C Using 
LQC/LTR Compensator 
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Figure 27. Singular Values of Closed-Loop Transfer 
Function G(Jco) K(ju) [I + G(Ju) K(Joj)H 
of Design C 



LOG FREQUENCY 

Figure 28. Stability Robustness Check for Design C 



Figure 28. Singular Values of the Sensitivity 
Transfer Function [I + GUto) 
of Design C 



Figure 29. Non-Linear Simulation of Response of MIMO System C to a Slow Temperature Reference Ramp 
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APPENDIX A 



APPENDIX B 


SISO SYSTEM A: 5 T = IN, N p Q M r NmrQth NtrJ 


AM.AZ.mjS 

* 3.64 

0.00 

0.00 

0.00 

0.00 

0.00 


B. matrix; 

46.20 



3-61 

- 1.09 

- 60.00 

0.00 

- 60.00 

0.00 


42.70 



0.60 

6.18 

10.00 

-6.20 

0.00 

0.00 


6.14 



0.00 

0.00 

9.06 

• 4.46 

0.00 

0.00 


0.00 



0.00 

6.29 

0.00 

0.00 

0.00 

- 6.29 


0.00 



0.00 

0.00 

0.00 

0.00 

186.00 

-.371 


0.00 



C MATRIX 
0.00 

1.00 

0.00 

0.00 

0.00 

0.00 





MIMO SYSTEM B: 

(SCALED) *T 

= (N g Np QmrNmrQtrNtr Vg| 





A MATRIX 
— 1.92 

0.00 

0.00 

0.00 

0.00 

0.00 

8.09 

B MATRIX 
11.47 

0.00 


3.24 

-. T27 

-.1217 

0.00 

-.012 

0.00 

- 13.6 

17.55 

0.00 


234.00 

2600 

-10.0 

-2600 

0.00 

0.00 

- 9.75 

1267 

0.00 


0.00 

0.00 

.018 

-.450 

0.00 

0.00 

040 

0.00 

0.00 


0.00 

31450 

0.00 

0.00 

0.00 

-31460 

0.00 

0.00 

0.00 


0.00 

0.00 

0.00 

0.00 

-.037 

- J7I 

0.00 

0 00 

0.00 


0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

- 10.0 

0.00 

10.0 


C MATRIX 
0.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 




1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 




MIMO SYSTEM C: 

(SCALED) 5 t 

= (N f Np QmrNmrQtrNtrT*.# 

v sl 




A MATRIX 
- 3.64 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

B MATRIX 
1639 1.01 

0.00 

7.84 

- 1.09 

- 1.212 

0.00 

-.150 

0.00 

0.00 

- 34.3 2.14 

0.00 

66.7 

259 . 

-10.0 

- 260 . 

0.00 

0.00 

0.00 

- 24.8 8J5 

0.00 

0.00 

0.00 

.1812 

-.460 

0.00 

0.00 

0.00 

0.00 0.00 

0.00 

0.00 

2516 . 

0.00 

0.00 

0.00 

-2616 

0.00 

0.00 0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

.465 

- J71 

0.00 

0.00 0.00 

0.00 

- 6.90 

0.00 

0.00 

0.00 

0.00 

0.00 

-240 

24.7 3.74 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

- 10.0 0.00 

10.0 

C MATRIX 
0.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 



0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

0.00 





